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Abstract 

We derive the dominant contribution to the large-distance decay laws of cor- 
relation functions towards their asymptotic limits for a spin chain model that ex- 
hibits both Haldane and Neel phases in its ground state phase diagram. The an- 
alytic results are obtained by means of an approximate mapping between a spin- 
1 anisotropic Hamiltonian onto a fermionic model of noninteracting Bogoliubov 
quasiparticles related in turn (via Jordan-Wigner transformation) to the XY spin- 
1/2 chain in a transverse field. This approach allows us to express the spin-1 string 
operators in terms of fermionic operators so that the dominant contribution to 
the string correlators at large distances can be computed using the technique of 
Toeplitz determinants. As expected, we find long-range string order both in the 
longitudinal and in the transverse channel in the Haldane phase, while in the Neel 
phase only the longitudinal order survives. In this way, the long-range string order 
can be explicitly related to the components of the magnetization of the XY model. 
Moreover, apart from the critical line, where the decay is algebraic, we find that 
in the gapped phases the decay is governed by an exponential tail multiplied by 
power-law factors. As regards the usual two points correlation functions, we show 
that the longitudinal one behaves in a "dual" fashion with respect to the transverse 
string correlator, namely both the asymptotic values and the decay laws exchange 
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when the transition line is crossed. For the transverse spin-spin correlator, we al- 
ways find a finite characteristic length which is an unexpected feature at the critical 
point. 

The results of this analysis prove some conjectures put forward in the past. 
We also comment briefly the entanglement features of the original system versus 
those of the effective model. The goodness of the approximation and the analytical 
predictions are checked versus density-matrix renormalization group calculations. 

PACS: 75.10.Pq Spin chain models; 03. 65.Vf Phases: geometric, dynamic or 
topological 

1 Introduction 

The Haldane phase |fl~), found in many low-dimensional spin systems, has attracted a 
great amount of attention in the last two decades both from the theoretical and from the 
experimental points of view. Its genuine quantum nature is signalled by two charac- 
teristic features. First, the excitation spectrum above the ground state (GS) displays a 
finite energy gap and, second, one can identify suitable long-ranged string correlation 
functions that measure a hidden topological order of the phase. The most intuitive idea 
to understand the physical features of the Haldane phase is probably the spin liquid pic- 
ture (21 : In a spin-1 chain with Heisenberg interactions and quantization axis directed 
along z, let us assign the presence of an effective spin- 1/2 particle with spin pointing 
up (down) if at the i-th lattice site Sf = +1 (—1) and no particles if Sf = 0. The 
Haldane phase is then interpreted as a liquid in which these effective particles carry no 
positional order along the chain but still retain antiferromagnetic (AFM) order in their 
effective spins. The positional disorder is associated with the absence of long-range 
order in the usual spin- 1 correlation functions 

C a (R) = {-1) R (S?S? +R ) , a = x,y,z 

whereas the spin- 1/2 magnetic order that we would get if all the sites with Sf = were 
taken off from the chain is measured by the asymptotic value of the string correlators 
@: 

O a {R) = (Sfe 1 ^^ s ? S? +R ) , a = x,y,z (1) 

for R — > oo. Interestingly enough, the Haldane gap has been interpreted as the exci- 
tation energy associated with a "spinon" (or kink) with respect to the hidden order (4). 
The nonvanishing values of the string-order parameters (SOP) 

O a = lim O a (R) 

R — 'oo 

can be understood as a spontaneous breaking of hidden (nonlocal) Z2 symmetries of 
the A — D Hamiltonian, as discussed thoroughly by Kennedy and Tasaki fT). From 
a numerical inspection of the string correlation functions (Q~|i computed on the first 
excited state with S* ot = 1, rather than on the GS, Elstner and Mikeska [4j argued that 
this excited wave function is characterized by a transition region with vanishing string 
correlations that connects two asymptotic limits with symmetry breaking and different 
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values of the hidden order. In a field-theoretic approach to spin- 1/2 Heisenberg chain 
1 5 1 such a kink is described as an effective particle - a soliton - moving with relativistic 
dispersion relation. When the system is moved away from criticality, due to the action 
of a relevant field, the soliton acquires a nonvanishing mass or an energy gap, in the 
condensed matter language. In ref. ]6] some of the authors have proposed a picture 
of the states that form the Haldane triplet at the isotropic point in terms of massive 
solitons and their bound states arising in the sine-Gordon formulation, valid in the 
neighbourhood of the critical line that marks the limit of the Haldane phase towards 
the large-Z) one (see below). The first solid numerical evidence of a nonzero Haldane 
gap has been provided by White and Huse |7J using the by now celebrated density- 
matrix renormalization group (DMRG) method. 

Actually, the Haldane phase is not restriced to spin-1 systems and can be found, 
for example, in spin-5 Heisenberg chains for every integer value of 5. According 
to ref. 10 the gap vanishes as the classical limit 5 — ► oo is approached as A a 
5 _1 exp(— irS) while the behaviour of the string order is more subtle: in order to have 
a nonzero value one has to generalize the string correlation function of equation (Q]i 
using not it in the exponential but 5-dependent optimal angles 6 n = (2n +1)tt/S with 
n = 0, 1, . . . , 5 — 1. Again, when 5 — > oo the resulting values of O a (9 n ) tend to zero. 

It is interesting to examine also how the features of the Haldane phase are destroyed 
by varying the parameters of the Hamiltonian out of the isotropic spin-5 Heisenberg 
model (5 integer). In this paper we shall stick from now on to the case 5=1 and 
consider two types of anisotropics along z: Ising-like interactions (parametrized by A) 
and single-ion terms (parametrized by D) 

H = Y^Si- S t+1 + (A - 1)5? 5f +1 + D(Sf ) 2 . (2) 

i 

The phase diagram of this model has been investigated in various papers with dif- 
ferent approaches J9l El [U . In order to fix the ideas we will refer to our recent deter- 
mination (6), reported (in a simplified form) in figure Q] Fixing a nonnegative value 
for A and varying D we encounter three gapped phases: the Large-Z? one in which 
O a = Va indicating the absence of magnetic order in the effective spin- 1/2 parti- 
cles. Their positional degrees of freedom are also disordered. In the Haldane phase 
the spatial disorder persists but magnetic order emerges. As a consequence both lon- 
gitudinal and transverse string order parameters (SOP) become nonzero: O a ^ 0. As 
pointed out in ref. iflOl on the basis of an exact solution for an integrable variant of 
(f2| with A = and in-plane anisotropy, the excitation gaps in the Large-D and in the 
Haldane phases have a rather different nature. Despite the fact that they are both found 
within the sector 5 t 2 ot = 1, the former corresponds to a flip of a single spin out of the 
xy plane while the latter is related to the breaking of a two-site singlet composing the 
resonant-valence -bond GS similar to the one of the spin-1 chain exactly solved by Af- 
fleck, Kennedy, Lieb and Tasaki [11]. Finally, by decreasing further the value of D, we 
pass in the Neel phase where both positional and magnetic degrees of freedom orders 
are signalled by a nonvanishing (spontaneous) magnetization along z 

M 2 Z ee lim C Z (R). 

ti — >OC 
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Table 1: Decay laws of string correlation functions defined in the text, according to 
Sect. HE of ref. [31]. To be compared with the results of this paper, including the 
explicit form of the algebraic prefactors, as reported in table [3] 



Phase 


O z {R) 


G H {R) 


Haldane 
Neel 


Expon. to ^ 
Expon. to 7^ 


Expon. to 
Expon. to 7^ 



At the same time O z ^ but O x , y = 0. Den Nijs and Rommelse (Q, Sect. HE) 
introduced yet another less-familiar string correlation function without spins at the ends 

and argued that Gh(oo) = in the Haldane phase but Gh(oo) 7^ in the Large-D 
and Neel ones. 

Hence we may select, equivalently, the pairs (O z , O x ) or (O z ,M z ) as order pa- 
rameters to classify the three types of behaviour. The universality classes associated 
with the two transition lines will be frequently denoted using the language of conformal 
field theory (CFT - see, for instance, |[T2l l9ll6lD. in particular by specifying the central 
charge c. We interpret the fully-disordered large-Z) phase with (O z = 0, O x = 0) and 
(O z = 0, M z = 0) as a spin gas. By crossing the c = 1 line we enter the Haldane phase 
where the effective spin- 1/2 experience a first magnetic ordering: (O z 7^ 0, O x 7^ 0) 
and (O z 7^ 0, M z = 0). Then, loosely speaking, at the c = 1/2 line the spin liquid 
crystallizes and the fully-ordered Neel phase can be interpreted as a spin solid with 
(O z ^ 0, O x = 0) and (O z ^ 0, M z ^ 0). Note the interchanged role of O x and M z 
(see below). In the Neel and Haldane phases Gh (oo) refers to the positional order of 
nonzero spins Q, so that it vanishes in the Haldane phase but G# (oo) 7^ in the Neel 
one. 

In order to determine the SOP numerically one has to extrapolate to the thermo- 
dynamic limit and to infinite distance the data computed on necessarily finite samples. 
However, apart from the qualitative statements made in ref. O about the exponential 
decay of the string correlation functions (as reported in table [TJ, the available litera- 
ture contains scarce information about the spatial behaviour of such correlators and the 
extrapolation may become problematic, especially close to the transition lines where 
the bulk correlation length becomes very large. In a particular case, namely the tran- 
sition from the Large-Z) to the Haldane phase, the low-energy physics is described by 
a compactified free boson field theory (c = 1 CFT). Once the compactification radius 
is known in some other independent way, one can read off the decay exponent of the 
string correlation functions from the set of scaling dimensions of the possible vertex 
operators. Interestingly, it turns out lfT3lfT4ll that, even if the lattice model has periodic 
boundary conditions (PBC), the vertex operators to be associated with string correlators 
belong to the sector with antiperiodic boundary conditions. 

The main purpose of this paper, instead, is to address the spatial behaviour of spin- 
spin and string correlation functions in the Haldane and Neel phases, making use of a 
solvable theory of spinless fermions. Starting from well inside the Neel phase, where 
the density of sites with Sf = is negligible, we approximate the problem by as- 
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Figure 1 : Ground-state phase diagram for the model (0 in the AFM region A > 0. The 
three phases are defined in the text. 



suming that the hidden magnetic order is frozen so that a given contribution to the GS 
wavefunction can be described by occupation numbers: no fermions if S? = and one 
fermion when \S*\ — 1, no matter the orientation, which is dictated by the underlying 
string order. The details of this approach will be presented in Sect. 12 actually it is very 
close to what done by Gomez-Santos in ref. lfl5l . The difference here is that we include 
also the single-ion anisotropy term and, in fact, the two formulations are related by a 
particle-hole transformation. The novelty is that we work out in detail the mapping of 
the spin-spin and string correlation functions (subsec. 12.11 ) onto fermionic correlators, 
so that we can derive in Sect. [3]the precise form of their asymptotic behaviour at large 
distances by exploiting the machinery of Toeplitz determinants. Sect. |4]reports a com- 
parison with DMRG simulations of the system in equation (f2]i while Sect. [5]is devoted 
to conclusions. 



2 Mapping onto spinless fermions 

The basic idea underlying the approximation used in this work is the spin solid picture 
of the Neel state(s): 

\N) = | TIT •■■ ITI •■■ U) (3) 

which is, in fact, the GS of the Hamiltonian (f2]l for A — > oo at fixed D or D — > — oo 
and A > 0. Actually the GS is doubly degenerate: for a given configuration of the type 
(0 with, say, Sf = 1 at the reference site i = 0, the energy is unchanged by the Z2 
transformation T = cxp(i7r J^j Sj) mat performs a 7r-rotation about the y-axis (spin- 
flip). We will refer to \N) and \N) = T\N) as Neel and anti-Neel states, respectively. 
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Now, in a perturbative fashion, when \D\ and/or A 3> 1 the effect of the transverse 
terms in the Hamiltonian Sf' y Sf+ 1 is to: 

i) create pairs of adjacent sites with Sf — Sf +1 = : | TJ.) — > |0 0); 

ii) move the zeroes in the AFM background, e.g.: | f 0} — ► |0 "f); 

iii) re-create a pair fj, or |f in place of a pair of adjacent zeroes. 

Notice that ii) preserves the AFM order, albeit not on nearest neighbours but mediated 
by string of zeroes (hidden order). Again, due to the AFM order (induced by A > 
and by the transverse terms), even if both states of iii) can be created in an "island" 
of zeroes, as far as the low-energy part of the spectrum is concerned, one of the two 
will be preferred according to the orientation of the surrounding spins, that is, by the 
hidden AFM order. Note also that \N) and \N) are connected through a large number 
of virtual processes, so that in the thermodynamic limit only one of the two will be 
selected by a spontaneous symmetry breaking mechanism induced by an infinitesimal 
staggered magnetic field. Alternatively, if the system under consideration is described 
by a thermal density matrix exp(— j3H), when (3 — > oo the GS reduces to a symmetric 
mixed state \N){N\ + \N)(N\. 

Once the question of the GS is accounted for, from the scenario above one can 
see that the orientation of spins with nonzero component along z is determined by the 
hidden order and can be taken for granted. The validity of such an approximation is 
ultimately measured by the values of the longitudinal SOP: the closer is O z to unity the 
higher is the AFM order of nonzero spins. We then introduce the following fermionic 
picture: assign a spinless fermion = q|— i) at site i if Sf ^ and no fermions 
|— i) = |0j) in the spin language if Sf = 0. (This notation for spinless fermions 
has a direct translation in the language of the XY model that will be introduced at the 
end of this section.) Process ii) is nothing but a hopping of spinless fermions, while 
processes i) and iii) represent annihilation and creation of pairs | + ; + The density 
of nonzero spins (Sf ) 2 is simply translated to the local fermion number rii = etc,, 
while, due to the underlying AFM order, the Ising-like term takes the form —AniUi+i 
that contributes with a negative energy when two fermions are present on adjacent sites. 

Hence, under the hypothesis of hidden AFM order, the dynamics of equation (O is 
reproduced by the following effective fermionic model 

Hi = X! (4 C J'+! + 4+i°o + c ] c ]+i + c :/+iCj - XrijUj+i + I) n,) (4) 

3 

(acting in a reduced Hilbert space Ji = <8>i7^ 2 ' 1 where Ti.^ denotes the local Hilbert 
space of a two-level system - as that of a spinless fermion or a spin- 1/2 introduced 
below). It should be observed that equation with D = is essentially equiva- 
lent (apart from an additive constant) to equation (2) of ref. |fT31 once a particle-hole 
transformation rii — > 1 — rii is performed at every site. 

Following Gomez-Santos lfl5l we now proceed to a further approximation on the 
fermionic Hamiltonian that is not amenable to an exact treatment due to the A-term. At 
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the Hartree-Fock level this term can be approximated as: 

n 3 n j+1 ~ (rij +n j+1 ) (n 3 ) - (ctcj+i(cj +1 c,-) + h.c.) + (ctcj +1 (c i+ iCj-) + h.c.J 

- (K) 2 - (4 +lC ,)(c]c J+ i) + (c, +lCj )( C ] c ] +1 )) 

where the expectation values (...) now are taken with respect to the GS of the quadratic 
Hamiltonian 

H HF = [C 1 + XA ) c h+i + C 1 - AB ) 44+1 + h-c. 
i 

+ (£> - 2An ) rij + A (rip - |A| 2 + |B| 2 ) (5) 

where the parameters 

n = (rij) , A = (ct +1 Cj) , S = (cj+ic,-) 

have to be determined self-consistently. The advantage of a Hamiltonian of the form 
(0 is that it can be diagonalized by means of a Bogoliubov transformation 

7] k = cos — c k + i sin — c'_ k 

where = c j ex P( — ij^) an d ^fe i s given by 

jg (cos fc — ft, + i7 sin fc) 

where 

h = 2A "" - g 7 =Iz** ( 6 ) 
2(1 + XA) ' 7 " 1 + AA W 



A k = \J (cos k — h) + 7 2 sin fc 2 

Note that, as we are interested in the thermodynamic limit, we do not specify here 
the boundary conditions on the spin and fermionic Hamiltonians. The momenta are 
quantized as Afc = 2ir/L and their precise location within the first Brillouin zone 
depend on the conditions imposed on the end sites. However, for L — > oo 



fc 



2tt 

dfc. 



Apart from additive terms of 0(L 1 ) the Hamiltonian in diagonal form is 



i?HF = 2(i + AA)^A fc Llm-\ 



u 
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and U = (D - 2\n Q )/2 + A(n„ - A 2 + B 2 ). In the thermodynamic limit, the self- 
consistent equations are 

1 1 f -Mnp, A) + cask 

The notation used in equation (|6]l is the one commonly used for the XY spin- 1/2 
model in a transverse field. In fact, by (inverse) Jordan- Wigner transform one gets lfl2l 



where the erf's are Pauli matrices at site j. This model is known to be critical at 
h = ±1 for 7^0, where it belongs to the c = 1/2 universality class (the same as the 
2D classical Ising model) and at 7 = for h € (—1,1) where the universality class 
becomes that of the compactified free boson, c = 1. 

From the numerical solutions of (IT])-® it turns out that in the Haldane and Neel 
phases A < and B < so that 7 > 1 (as long as \\A\ < 1 - some representative 
cases are listed in tabled, while most studies are limited to |7| < 1. As a consequence 
the region j 2 < 1 — h 2 corresponding to oscillations with wavenumber different from 
7r lfl6l FPU is not present in our case. However, having 7 > 1 does not affect the 
critical condition we are interested in, that remains h = ±1. In these cases we have 
just c = 1/2, as reported before |9| for the Haldane-to-Neel transition. At A = 0, 
D = — 2 this transition line merges with the boundary towards the so-called XY phases 
corresponding to c = 1, Interestingly this change of universality class is captured 
also by our approximation since for A = 0, D = — 2 the self-consistent solution 
yields just 7 = 1 and h = 1. From the data in table |2] one can also estimate, for 
example, the critical value of D at fixed A = 1; the result is D c = —0.214, which 
is not in good quantitative agreement with the numerical value D c = —0.315 ll9l [T4l . 
The perturbation of the isotropic Heisenberg Hamiltonian with A > 1 and D = 0, 
instead, seems to be better described by the spinless fermions approach; already at 
this level of approximation the value A c = 1.125 found in |fT31 is close to our best 
DMRG independent estimate A c = 1.1856 ifTTl . Even if it is likely that the inclusion 
of configurations with nearest-neighbour parallel spins could improve the results, as 
discussed by Gomez-Santos lfT31 . we do not insist along this line here because we 
are ultimately interested in the decay laws of correlation functions that are essentially 
dictated by the universality classes. In fact, it is important to stress that neither the 
extension to D ^ 0, nor the extension of the model as in equation (8) of ref. |[T5l 
modify the universality class of the transition, that remains of the c = 1/2 (or Ising) 
type for A > 0. Although the location of the critical points and of the prefactors 
depend on the values of the parameters, the scaling dimensions of the operators in the 



Table 2: Self-consistent estimates of the three decoupling parameters no, A and B of 
equations ©-(O for some choices of A and D in the Haldane and Neel phases. It must 
be kept in mind that the continuum versions of the self-consistent equations neglect 
some 0(I/ _1 ) terms coming from isolated contributions at wavenumber or 7r. Last 
two column contain the corresponding parameters h and 7 of the effective XY model 
according to equation ©. 



A 


D 


n 


A 




B 




h 




1 




1 





0.709 


-0 


.158 


-0 


.253 





.842 


1 


.49 


1 


-0.125 


0.745 


-0 


.137 


-0 


.246 





.936 


1 


.44 


1 


-0.200 


0.774 


-0 


.117 


-0 


,240 





.990 


1 


.41 


1 


-0.250 


0.800 


-0 


.0979 


-0 


.235 


1 


.03 


1 


.37 


1 


-0.300 


0.816 


-0 


.0866 


-0 


.231 


1 


.06 


1 


.35 


1 


-0.315 


0.820 


-0 


.0837 


-0 


.230 


1 


.07 


1 


.34 


1 


-0.330 


0.824 


-0 


.0811 


-0 


.229 


1 


.08 


1 


.34 


1 


-0.345 


0.828 


-0 


.0786 


-0 


.228 


1 


.09 


1 


.33 


1 


-0.400 


0.841 


-0 


.0706 


-0 


.223 


1 


.12 


1 


.32 


1 


-0.450 


0.850 


-0 


.0645 


-0 


.219 


1 


.15 


1 


.30 


1 


-0.750 


0.893 


-0 


.0406 


-0 


.198 


1 


.32 


1 


.25 


1 


-0.875 


0.904 


-0 


.0344 


-0 


.190 


1 


.39 


1 


.23 


1 


-10 


0.996 


-0 


.000317 


-0 


.0433 


6 


.00 


1 


.04 


■3 


-0.125 


0.991 


-0 


.000853 


-0 


.0649 


5 


.04 


1 


.33 



continuum field theory (i.e. the decay exponents of the correlation functions) do not 
change when we move along the c = 1/2 line. Nonetheless, due to the lack of an 
explicit mapping of the spin-1 strings onto the corresponding correlators in the Ising 
fermionic field theory, up to now the exponents appearing in the large-distance decay of 
string correlation functions were unknown. This is precisely the subject of subsecs. 13.11 
and |3.3l Eventually, we note that alternative pictures of the Haldane gap in fermionic 
language can be derived by perturbation theory near the Babujian-Takhtajan integrable 
biquadratic spin- 1 chain |[T8l or from two-leg ladders with ferromagnetic coupling on 
the rungs |[T9l . 

At this stage it is interesting to compare the entanglement properties of the original 
spin-1 model (eq. (0) with those of the XY spin- 1/2 chain resulting from the mapping. 
On the one hand, for the former it has been shown 120 1 that at the isotropic Heisenberg 
point A = 1, D = there is long-distance spin-1 (qutrit) entanglement in the ther- 
modynamic limit for two sites arbitrarily far apart. It is reasonable to expect that this 
entanglement survives in a neighbourhood of the isotropic point. On the other hand, 
in ref. ED it is stated that the qubit entanglement in the XY model with transverse 
field vanishes beyond a distance of order 7 -1 . In our case 7 > 1 and the degrees of 
freedom of the qubits represent the presence or the absence of an effective particle with 
I I = 1. Therefore we are led to speculate that wherever there is full spin-1 entangle- 
ment in the vicinity of the Heisenberg point, this is due to the spin correlations between 
the sites with Sf 7^ 0. Recalling the hypothesis of underlying string order and imag- 
ining to eliminate the sites with Sf = 0, the qualitative picture of the long-distance 



9 



entangled states in the Haldane region if that of a Greenberger-Horne-Zeilinger state 
flSl with effective AFM order | . . . Titl ...) + !•■■ IT4T • ■ ■ > ■ 



2.1 Mapping for the spin-spin and string correlators 

We shall exploit now the mapping from spin-1 to spinless fermions, based on the exis- 
tence of an underlying string order, to translate the various spin-1 correlation functions 
onto expectation values of strings of fermionic operators that can be computed exactly 
when the Hamiltonian has the form ©. Let us start from the z-component of the spin 

- mo - — ^-tfj (id 

where Kj = exp(i7r^ i<; ■ ni) = Yli<j(~ a i) * s a Jordan-Wigner tail that accounts 
for the correct sign when S| 7^ assuming, conventionally, that the first nonzero spin 

is pointing up. By inserting the expression S*| = l+ ° 3 Yli<j(^' J j) mt0 tne defini- 
tion of the longitudinal spin-spin correlation function and using the properties of Pauli 
matrices one finds 

CAR) - ±(-l) R ((l + a])H(-ai) J] (-a* k )(l + a] +R )) 

k<j k<j+R 
I j+R j+R j+R-1 j+R-1 \ 

= 4 <n a ^ + < n o + < n + < n )■ ^ 

\ k=j k=j + l k=j k=j+l J 

As far as the transverse correlation functions are concerned, it can be checked by 
direct inspection on a generic configuration with perfect string order that the identifi- 
cation 

j+R-1 , z \ 

CM ^( VV» • S I S I+r) •< II (-2^)^ (I 3 ) 

k=j+l ^ ' 

has the correct action, since the only cases in which the l.h.s. does not break the string 
order are those with S k = 0, that is <r| = —1, on all sites between j and j + R. 
The product on the r.h.s. of ( fT3T > is exactly the expression involved in the so-called 
emptiness formation probability (see, for example, ll23l and refs. therein). 
Let us now study the spin-1 strings. Along the z-direction we have simply 

e -5: &<3 SJ = JJ (1 _ 2 {Slf) - J] (-°k) • (14) 

k<j k<j 

Again by using the relation 5J = Y[i<j(~ a j ) an d plugging me string written 
above into eq. (Q]i one gets 

o.( J B)-/(^)n(-^ ,+ ff 1 (-^) n ^C-^y 13 -)) (i5 > 

\ ^ ' k<j k=j+l k<j+R ^ 'I 



; (^) (^ ± ^)> = -i(i + K> + K + «) + (^l + «)) d 6 ) 
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Note that in the language of the effective XY model, the Neel correlation function 
( fT2b involves a string of Pauli operators whereas the string correlation function ( TToT l 
involves only one- and two-points correlators of the er's. Thanks to equation ( TPfl i we 
easily obtain also the pure-string correlation function as: 



li+R 



G H (R)^(-l) R+1 {J[a]j. (17) 

From equation ( fT2l ) we see that, in this approach, Gh (R) is nothing but the first term 
of the usual spin-spin correlation function C Z (R) apart from the prefactor. 
Along the x-direction, instead, we exploit the fact that 




that is, apart from an overall sign, the operator above performs a swap between Sj = 
1 and S? = — 1 leaving the Sj = component isolated. This swap is important 
because it can be checked by direct inspection that, for every possible combination of 
the spin at sites j and j + R that respects the hidden AFM order, both with an even 
or an odd number of nonzero spins in between, the action of O x (R) produces only 
one allowed configuration and some other forbidden ones. More precisely, using the 
spin- 1/2 operator a x that changes empty sites into occupied sites and viceversa we can 
write 

O x {R) - (-l)"- 1 ^ 

where the inner spin-1 transverse string contributes with the sign prefactors. Thanks 
to hidden order in our reduced Hilbert space, the spin- 1/2 configurations generated by 
a x represent the allowed spin- 1 states and the forbidden ones are automatically filtered 
out. The coefficient 1 /2 comes from the matrix elements of S x at sites j and j + R. 

Now, thanks to the fact that the Hamiltonian © is quadratic in the fermionic oper- 
ators, all the correlation functions can be evaluated using Wick's theorem. Following 
the notation of the seminal paper by Lieb, Schultz and Mattis [24] we introduce the 
operators Aj = c ] + c j an d Bj = cj — Cj that allow to express the basic two-point 
correlations as 

(°f°m) — (BiA[ + iBi + i ■ ■ ■ A m -iB m -iA m ) 

(of O = {AiBiA m B m ) 

with Q lm = (AiA m ) = Si m and Si m = (BiB m ) = -Si m . If we further assume 
translational invariance (i.e. PBC) along the chain we have 



\ i o TO 



G-i G-2 ■ ■ ■ G, 



l — m 



Gm-l-2 ' ' ■ G_i 



(18) 



Gq — G m -iGi- m (19) 
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where G-r 

2(n 



= (B 3 A 3 _ 



= <(cJ 



d(4 



r) = -(Aj +R Bj). In particular, G - u --j ~ ^A-, 
cj)) = 2{rij) - 1 = (<r|), independent of j and (VjVj +R ) = (o|) 2 - GrG-r. 
The independence of O z (i?) and O x (R) is given directly by (<t|(t| + jj) and (crja'j_ ( _ fi ) 
respectively. The ordinary correlators C XtZ (R) require a step more since they involve 
strings of Pauli operators. For example, each of the terms in equations fl% and (JT7J 
has the form (Y[ k B^Ak). When R — > oo all the four terms in equation (fT2l tend to 
coincide so that 



C Z (R) — (— 1) Gh{R) — (BjAjBj + iAj + i ■ ■ ■ Bj + R-iAj + R-\Bj + RAj + R) 



Exploiting Wick's theorem, Caianiello and Fubini [25] have shown that the expectation 
value above can be expressed as a Pfaffian 



Pf|5_i S- 2 



S-R+l 
S-R+2 



S-R 
S-R+l 



Go 



G-R+l 
G-R+2 



G-R 
G-R+l 



S-! 


S-2 


Gr-2 


Gr-3 ■ 


■ G_i 


G-2 






Gr-x 


Gr-2 ■ 


Go 


G-x 






Gr 


Gr-x ■ 


Gx 


Go 








Q-x ■ 


Q-R+X 


Q-R 



Q- 



Q-2 

Q-i 



Thanks to the fact that Qi^ r , 
nant ll26l and we get 



Si-t m — this Pfaffian reduces to a Toeplitz determi- 



C Z (R) = 



—Go 
-Gx 

-Gr-x 

— Gr 



-G-x 
—Go 

— Gr-2 

—Gr-x 



-G-r+x 

-G-R+2 

—Go 
-Gx 



—G-r, 
—G-r + x 

-G-x 
—Go 



(20) 



So, the determinants of the matrices with entries Gj becomes the central quantities of 
our analysis. 

The matter is more complicated for the transverse spin-spin correlator essentially 
because it eventually involves a Toeplitz determinant generated by a matrix-valued 
symbol that may also become singular. According to ref. l27| this case is not yet 
solved in the theory of Toeplitz determinants and in ref. l28l it has been suggested 
to extend directly the procedure valid in the nonsingular case. Fortunately in our case 
a workaround is possible: thanks to a suitable diagonalization, we are able to com- 
plete the calculation of the dominant contribution to C x (R) in terms of a product of 
Toeplitz determinants, each one computed using the Fisher-Hartwig conjecture (see, 
for instance, App. A. 2 in 11231 ). The details of this procedure are reported in the Ap- 
pendix. 
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3 Asymptotic decay laws 

3.1 Longitudinal string correlation O z (R) 

The first object we will compute is the longitudinal string correlator. From equations 
( fT6T > and ( [19] l we get 

O z {R) = -\ [(l + (a~)) 2 -G fl G_ 

Following Barouch and McCoy |fl6l (in particular their eq. (6.12)) we express Gr as 
follows 



1 



2tt 



2jt 



G R = - — / dke~ ik ( R+ V 



(1 - \- x e ii: ) (1 - 1 e ifc ) 
(1 - aYV^) (1 - A^" 1 e _ifc ) 



-,1/2 



with 



and 



c(c Ifc ) = c- 



h± x /h 2 -(l-j 2 ) 
1^7 ' 



1 f 2 " 
— / dke 

(21) 
(22) 



1 (1 - Aj-y fc )(l - Aa V fc ) 

(i - Aj-y- ife )(i - A^y- ife )' 



Note that since 7 > 1 the two roots of the numerator are always real; the behaviour for 
R — ► 00 is controlled by A2. From equations (6.17), (6.14) and (6.18) in ref. Ifl6l we 
have, respectively: 



• Haldane phase h < 1(A 2 > 1) 

O z (R) ~ O z 

• Critical line h = 1 (A2 = 1) 



1 e~ 2R /Z 
8^ R 2 



, £=l/mA 2 



• Neel phase h > 1 (0 < A 2 < 1) 

1 e" 2i? /« 
+ - e=-l/lnA a 

07T ft-* 

In every case the asymptotic value O z 7^ is simply interpreted as a non-saturated 
value of the magnetization along z in the XY model in transverse field 



O z = - 



(23) 



where (erj) = Go(h, 7) can be computed using equation (f2Tb at R = 0. 
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3.2 Longitudinal spin-spin correlation function C z ( R) and pure string 
correlator Gh(R) 

The asymptotic behaviour of the Toeplitz determinant in equation (l20t can be found 
using the same technique as in 0291 , since (apart from a sign) the generating symbol 
c(e lfc ) is essentially the same used by Wu. Then we find: 

• Haldane phase h < 1 (A2 > 1) 

C Z (R) * (-1) R ^G H (R) = -L(l-Ar 2 ) 1 / 4 (l-A 2 - 2 )- 1 /4 (1 _ Ar i A2) -i/2^ 

which corresponds to the known decay behaviour at the isotropic Heisenberg 
point, as predicted by the nonlinear er-model approach (see, for example, l30l ). 
Moreover, in refs. OTl [32l it was argued that the same behaviour of the con- 
nected longitudinal correlation function persists also in presence of a staggered 
magnetic field; in this sense such a behaviour could be considered a signature of 
the Haldane phase, robust against anisotropic perturbations. 

• Critical line h — 1 (A2 = 1) 

C Z (R) ~ (-l) R+1 G H (R) = e^Vi^-a^^ 

where A = 1.282427130 . . . denotes Glaisher's constant 11291 . 

• Neel phase h > 1 (0 < A 2 < 1) 

C Z (R) ~ (-1) R+1 G H (R) = (l-A r 2 ) 1 /4(i_ A 2 ) i/4 (1 _ Ar i A2) -i/2 



-2R/V 



2^(A 2 T 1 - A 2 ) 2 R 2 



Apart from the nonzero asymptotic value for h > 1, that serves as an order parameter 
for the Neel phase (the ordered phase T < T c in Wu's paper l29l ). it must be noticed 
that both the power of R in the denominator and the exponential constant are different 
on the two sides of the transition. The roots A12 and the bulk correlation length £ are 
the same as in subsec. I3.1l (see eq. d22l). 



3.3 Transverse string correlation function O x (R) 

• Haldane phase h < 1 (A2 > 1). The nonzero asymptotic value O x comes 
from the long-range order lim/^oo (ctqCt^) in the XY model with spontaneous 
breaking of the symmetry a x — > — a x . The result can be borrowed directly from 
equation (4.1) of lfl6l 



O x (R) 



[7 2 (l-fr 2 )] 
1 + 7 



1/4 



1 



-2R/e 



2ttR 2 (A a - A^ 1 ) 2 



(24) 



with £ having the same meaning of subsec. 13.11 
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• Critical line h = 1 (A2 = 1). There is no long-range-order in ((JqCTr), that 
decays to zero as i? -1 / 4 as expected from the scaling dimension 1/8 of the 
primary operator in the c = 1/2 CFT lTP2l . Using equation (4.7) in IfTBI we 
have 

O x{ R) - -7^^ 2 V^-3_i_ (25) 

• Neel phase h > 1 (0 < A2 < 1). Equation (4.25) in ref. Ifl6l is translated to 

O x {R) ~ -^=^T [(1 - AD-^l - A r 2 )(l - K 1 ^ 1 ?] 1/4 ■ (26) 

We should stress that the critical exponent in equation (l25T l differs from the one 
in equations (124-b and (|26*1 i: it is not possible to recover the decay behaviour at 
h = 1 from the functions found for h > 1 or h < 1 simply by letting — > 
in the exponentials. Qualitatively, the reason is that the correlation functions 
should be described by a unique scaling function T(r) of the variable r = i?/£, 
but the asymptotic expansions in the off-critical regime and in the critical regime 
are different. The former corresponds to r ^> 1 while the latter to r — > for 
any large but finite value of R. A similar argument holds also for the longitu- 
dinal spin-spin correlation function C z (R) of the previous subsection. Although 
possible in principle (see, for instance, l33l for the 2D classical Ising model), 
the derivation and the usage of the whole scaling functions is beyond the scope 
of this paper. Finally, we notice from the equations above, as compared to those 
of subsec. ( 13.2b . that the correlators O x (R) and C z (R) play a dual role above 
and below the transition line; when one order parameter is vanishing, the other 
is not. Here we do not have an explicitly duality relation between order and dis- 
order lattice operators for the spin-1 model as in the Ising case (see, however, 
ref. l34l for the XY chain). Hence, what is a nontrivial fact is to see that also the 
decay laws interchange when the transition line is crossed. 



3.4 Transverse spin-spin correlation C x (R) 

From the analysis reported in the Appendix, we can prove a conjecture already put 
forth in ref. fl"5l . namely that the transverse correlation function decays always ex- 
ponentially even when one crosses the critical line. Here we can be more precise and 
derive also the power-law terms in front of the exponential 

cxp(-i?/S) ^_ 2 
C X (R) — , ~=J^^ (27) 

where (3 and (3' in the Haldane phase, along the critical line and in the Neel phase 
take, respectively, the form written in equations d37| i, (l4TT i. d38l ), ( l42l . ( T39b and d43l 
in the Appendix. In particular we have checked that both (3 C and (3' c for h — 1 and 
7 > 1 are nonzero. Hence, despite the fact that the system is critical, the transverse 
correlation function exhibits a. finite characteristic length H. As far as the exponent r\ x 
is concerned: 



15 



• Haldane phase h < 1 (A2 > 1) and Neel phaseh > 1 (A2 < 1): T] x = 1/2; 

• Critical line h = 1 (A2 = 1): rj x = 1/4. Despite the fact that 5(/i = 1) < 00, 
the algebraic prefactor is the same power-law that describes critical correlations 
in the quantum Ising model. 

From the values of h and 7 reported in the last two columns of table [2] we have com- 
puted S[h(X, D),j(X, D)]; for example when A = 1 we find that 3 decreases steadily 
as D is decreased towards larger negative values, passing from the Haldane to the Neel 
phase. This behaviour is consistent with the numerical best-fit estimates of S made in 
the next section. 



4 Comparison with DMRG results 

The results of the previous section regarding the long-distance decay of ordinary and 
string correlation functions are summarized in table [3] where: 

fo(R)=A ^t^l (28) 

M R)=A 2+ A ^t^M (29) 
V R 



f 2 (R)=A 2 + A ^t^m . (30 ) 

Within the approximation of hidden order and for large R these asymptotic laws are to 
be considered exact and valid for the Haldane and Neel phases and associated transition 
line as specified in table [3] It should be noted that fi(R) and f 2 {R) agree with the 
general form argued for the <i-dimensional Ising model (see, for example, 11351 and 
refs. therein) although the derivation of the latter did not include the case of string 
correlation functions. 

The two functional forms fx ,2 now can be used to extract, for example, the asymp- 
totic value of string order correlation functions computed numerically; in this sense 
Aq, Ax and A2 may be regarded as free fitting parameters. The goodness of the best-fit 
procedure can be assessed by computing the reduced x 2 '- 

o Y] squares of differences 

X = . 

# of data points — # of fit parameters — 1 

Clearly one could use many other different functions to extrapolate the correlators 
to R — > 00. However, as recalled in the introduction, the literature contains very few, 
empirical, information about the asymptotic approach to the limit values of the string 
correlation functions. Our study was motivated by this fact and so here we perform 
a comparison between /o, fx and f 2 by examining their capability to fit the spin-spin 
and string correlations evaluated numerically through the DMRG. Actually, following 
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Table 3: Expected asymptotic behaviour of string (O) and usual (C) correlation func- 
tions in the Haldane and Neel phases of model (O and along the critical transition line 
separating them. The fitting functions /o,i,2 are defined in equations (l28l>- (13011. Note 
the interchanged role of O x (R) and C Z (R) above and below the transition line. 



Phase 


C.f. 


Decay law 


Haldane 


c z 


h(A 2 = 0) 


Transition 


c z 


fo (Ai 1 = 0) 


Neel 


c~ 


h 


Haldane 


c. 


h(A 2 = 0) 


Transition 


c 


fo 


Neel 




fi(A 2 = 0) 


Haldane 


o z 


h 


Transition 


o z 


h (V = 0) 


Neel 


o z 


h 


Haldane 


o x 


h 


Transition 


o. 


fo (^r 1 = o) 


Neel 


o x 


fi(A 2 = 0) 



the idea of ref. j30l . in order to take into account the PBC we employ the left-right 
symmetrized expressions of equations (f28T>-(f30b 

F t (R)= MR) + f * L - R) , t = 0,1,2, (31) 

at least for the correlation functions in the ^-channel. As regards O x , while transla- 
tional invariance implies that it depends on the difference between the sites at the ends 
of the string, it is not always guaranteed that it depends only on the distance on the 
ring. In other terms the expectation value 

(5fe i7r ^=*+i s *5J) 

may differ from the same expression with i and j interchanged. In fact, using the 
properties of the exponentials of spin-1 operators, exp(mSf) = cxp(—inSf) and 
Sf exp(mSf) = —Sf, it can be shown that the expression above can rewritten as 

where S* ot = J2i=i The point is that in general the GS of an anisotropic spin chain 
is not invariant under the action of cxp(i7r5'f ot ) and so a direct inspection is required 
case by case in order to decide if a symmetrized fitting function has to be used or not. 

The asymptotic limits (i.e. the values of A 2 ) resulting from of a series of best- 
fits made on DMRG data obtained by fixing A = 1 and letting D to vary across the 
Haldane-Neel transition from —0.125 to —0.875 are plotted in figure [2] It is seen that 
the nonvanishing values of O x and M z characterize, respectively, the Haldane and the 
Neel phase. It is reasonable to expect that the location of the critical point as the value 
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D 



Figure 2: Asymptotic values (order parameters) attained by C z (dots), O z (up triangles) 
and O x (down triangles). The limits to R — ► oo correspond to the values of the best-fit 
parameters A2 for the fitting function that, case by case, gives the smallest value of 
X 2 . The empty triangles result from fitting the transverse string correlation functions 
on half chain (see text for explanation). 

of D at which the two order parameters vanish leads to two slightly different estimates. 
However, with more accurate methods the critical point was previously found to be 

d c = -0.315 mini. 

In the simulations we have fixed the total length of the chain to be L = 100 sites 
and computed the GS properties by retaining from 243 to 324 DMRG states in the 
sector with S^ ot = 0, which is the only good quantum number that we could use. 
All the functional forms derived above are asymptotic so we cannot expect them to be 
reliable for very short distances. Therefore, we have conventionally excluded the data 
with R < 5 from the fitted points. In the Neel phase the GS tends to become doubly 
degenerate in the limit L — ► 00; in order to take into account this difficulty we have 
built the reduced density matrix by targeting the two low-lying states rather than just the 
GS. Finally we have performed three finite-system sweeps to achieve a better accuracy. 
In the cases we have considered, the transverse string correlation O x (R) turned out 
to be symmetric with respect to the middle of the chain except for D = —0.75 and 
D = —0.875. For this reason we have repeated the fit using directly the functions 
of equations d28t . d29l ) and d30b without symmetrization selecting only the points in 
the first half of the chain. The asymptotic values are essentially unaffected, with the 
exception of those referring to the critical point. In general when the results of the fit 
are such that A\ ^> L (typically close to criticality) we conclude that the exponential 
tail of the fitting function is essentially saturated to unity and an algebraic fit would 
produce the same result. 

As far as the best-fitting functions for C z (R) and O x (R) are concerned, the passage 
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from the type f\ to the type [2 going through fo at the critical point, as in table [3j 
actually takes place gradually in the interval D E (—0.345, —0.315), the worst values 
of x 2 being of order 10 -5 . The best choice to fit the transverse spin-spin correlation 
function C X (R), instead, follows the prediction of table^^i except at the critical line, 
where it becomes Fq) with a deviation \ 2 < 10 -8 . Finally, the longitudinal string 
correlator O z (R) is very well fitted by F2, in agreement with table[3] with \ 2 ~ 10 -9 
or better. 

It is also important to check quantitatively the goodness of the Hartree-Fock ap- 
proximation. The decoupling parameter in the fermionic version are no = (nj = o), 
A = (c\cq) and B = (ciCq) where we have selected a reference site "0" invoking 
translational invariance. In the original spin-1 formalism it can be checked directly that 

n Q = ((S~) 2 ) , A = 1 -{SI (5 +5r + S+So) S* ) (32) 

The operator ciCq destroys a couple of fermions in adjacent sites; in the spin lan- 
guage they could be fj. or J,| depending on the surrounding sites in order to respect 
the AFM order. Let us express the GS in the form |GS) = a\ ft) + f3\ JJ.), where | ft) 
denotes a linear combination of states in which the first nonzero spin along-z is di- 
rected upward and | JJ,) the same state with all the spin reversed. Only one of the terms 
in (So Si + S^So) will act on | ft) respecting the AFM order and the other term will 
thereby act on | JJ.). When the scalar product with (GS| is taken, the states from | ft) 
will not mix with those from | JJ.) . Therefore we try with the expression 

B = -\{{SoS+ + 3+Sr) S$Sf). (33) 

In table |4] we report the values of the decoupling parameters for a set of points in 
the Haldane and Neel phases, comparing the DMRG values with the numerical so- 
lution of the self-consistent equations using 100 iterations from different choices of 
initial conditions. Having the DMRG estimates for no and A we may also produce 
a "hybrid" estimate of the critical point by setting ft, = 1 in equation (|6]l and then 
solving for D C (X) = 2[A(n DMR G — A DMRG ) — !]■ With A = 1 we find for ex- 
ample D c = —0.254, that compares slightly better than the fully-self-consistent value 
(D c = —0.214) to the accepted numerical one D c = —0.315. 

Apart from the value no, which quantifies the number of spins with nonzero pro- 
jection along z, we expect that the goodness of the mapping used in this work is higher 
when the hidden order is larger. Therefore, as a final check, we have repeated the pas- 
sages of Sect. |2](see |[36l for details) including also a biquadratic term i *^i+i) 2 
in the spin-1 Hamiltonian (0. For A = 1 and D = the ground state of the model can 
be found exactly ifTTI fTl using the valence-bond picture: each spin-1 is viewed as the 
triplet sector of a pair of spin- 1/2 particles and the ground state is constructed by cre- 
ating a sequence of singlets between adjacent sites. In this case the string correlation 
functions can be computed exactly and it turns out that O x ^ z (R) = —4/9 independent 
of R. At the isotropic point with biquadratic term the self-consistent equations are 
solved by no = 2/3, A = B = —2/9 and the effective parameters of the XY model 
become h = 3/5 and 7 = 4/5. From equation (I22I 1 we find Ai = A2 = 3 so that in 
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Table 4: DMRG (L = 100) versus self-consistent (s-c) estimates of the three decou- 
pling parameters uq, A and B of equations (l32l and ( f33l >. It must be kept in mind that 
the continuum versions of the self-consistent equations neglect some 0(L~ l ) terms 
coming from isolated contributions at wavenumber or x. 



A 


D 


no.DMRG 


^Ldmrg 


-Bdmrg 


«0,s-c 


c 


B s -c 


1 





0.667 


-0.166 


-0.30080 ±0.00005 


0.709 


-0.158 


-0.253 


1 


-0.125 


0.702 


-0.151 


-0.28908 ±0.00005 


0.745 


-0.137 


-0.246 


1 


-10 


0.996 


-0.000324 


-0.0442 


0.996 


-0.000317 


-0.0433 


5 


-0.125 


0.991 


-0.000860 


-0.0654 


0.991 


-0.000853 


-0.0649 



equation (l23l we have just (<r|) = 1/3 and O z = —4/9. Interestingly enough, even 
if the XY model does not have an explicit rotational symmetry as the original spin-1 
Hamiltonian, by inserting these values of h and 7 into the constant part of equation 
( 1241 we find again O x = —4/9. This accordance can be taken as a positive check of 
our approach. 

5 Conclusions 

In this paper we have reconsidered and extended the approach of ref. |fl~5l to the GS 
properties of spin-1 anisotropic quantum chains. We have included a single-ion term in 
the Hamiltonian and, moreover, we have analyzed explicitly how the spin-1 correlation 
functions are written in the spinless fermions language and then in the framework of the 
XY model in a transverse field for effective spin- 1/2 degrees of freedom. In particular, 
we have focused on the decay laws of the string correlators towards their asymptotic 
values which apparently were missing in the literature. 

The decay laws of string and spin-spin correlation functions (in the longitudinal 
channel) are all related to the generating function c(e lfc ) of equation ( f2ll and the deter- 
minants of the Toeplitz matrices derived from it. The asymptotic behaviour of the trans- 
verse correlation function C x (R), instead, originates from a product of two Toeplitz de- 
terminants (see the Appendix, in particular eq. d35ll). The leading terms in the regime 
R ^> 1 for the various correlators are discussed in Sect. [3] and summarized in table [3] 
In brief the most interesting points unveiled by the approach used here are: 

• The nonvanishing string-order parameters of the spin-1 model (ffj) are simply 
interpreted as the magnetization along x and z in the XY chain with transverse 
field (eqs. <|23]) and (EH). 

• There exists dual behaviour of O x (R) and C z (R) above and below the transition, 
both for the asymptotic order parameters and for the decay functional forms. 

• The explicit calculation of C x (R) allows us to prove an unusual feature in statis- 
tical mechanics, already conjectured by Gomez-Santos 1131 : the spin-spin trans- 
verse correlation function exhibits always a finite characteristic length 3 (eq. 
d27l i) even when the system becomes critical. 
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The analytical results are supported by comparison with a numerical (DMRG) study of 
the model, especially for the correlations C x and O z . A more detailed comparison be- 
tween the analytical and the numerical estimates should take into account: i) finite-size 
effects due to a finite total length L while in Sect. [2]we passed readily to the thermody- 
namic limit; ii) corrections for finite distance R beyond the dominant ones. Although 
in principle they can be computed systematically, in this paper we have limited our- 
selves to the leading terms in order to derive analytical expressions with the smallest 
possible number of fitting parameters. 
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Appendix: Toeplitz formulation of C X (R) 

The fermionic version of equation dot reads 

j+R-l 



j+R-l 

-{S+S- +R + SjS+ +R ) = {A J \{{l-2n k ) Yl (l-n fc ) J] (l-2n fe )^ +fl ) 

k<j k=j+l k<j+R 

'j+R-l \ I j+R-l 

(Bj\ n c k ci\ A j+R ) = (4 1 n c k4\c] +R ) 

,k=j+l J \k=J+l 

( j+R-l \ I J+R-l \ / j+R-l \ 

+^ c ] n ck °i c 3+n) - fa n ° kc i c \+r) - (°i n ck °i ^+r) 

\k=J + l J \k=j+l J \k=j+l ) 

By observing that 

' j+R-l \ /( I j+R-l \ \\ ( j+R-l 

<s f ( n c 4 H+* )= ( u n ) c I+r) )=-<<* n 

K k=j+l J \\ \k=j + l J J I \k=j + l 

' ' j+R-l \ // (j+R-l \ \ f \ /j+R-l 

{c H n ckc i\ c J+R) = (\ c ] [ n ck °i ) c j+r\ )=-( c n n c *4H+i 

k=j+l j \\ \k=j + l J ) I \k=J + l 

we can write C X (R) = {-l) R {SfSj +R + SjS+ +R )/A (using translational and U(l) 
rotational invariance about z) as 



C X (R) = --v/detMi - VdctM 2 (34) 
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where the two terms come, respectively, from the Pfaffians 



Pf|iFi iF 2 
iFi 



iF R - 3 \F R -2 



iF 1 



iF 2 



-H-i 

-Ho 

-Hr-a 

-Hr-3 
-Hr-2 



-H_ 2 
-H-i 

— Hr-5 

— Hr-4 

— Hr-3 

-iFi 



—H_r + i —H-r 
—H-r+2 — H-R+l 



-#-2 

-H_ x 

—\Fr_2 



and 

Pf|iFi LF a 



iFr-2 
iF 2 



-H_ x 
-Ho 

-Hr-4 
-Hrs 
-Hr-2 



-H_ 2 

-H-i 

— Hr- 5 

— Hr-4 
-Hr^ 3 

-iFi 



-H-2 

-H-i 

-Ho 

-iFr-2 

-iFi 



with 



F, 



i-j = WCl) 



2^ Jo 



dke-W-j) f ( e ife ) 



1 /" 

7 sin A; 



e ife ) e" ife 



2 J (cos A; - h) 2 + i 2 sin fc 



, h(e ik ) 



-H- 3 

-H-2 

-H-i 
-iF R -i 

-iF 2 
-iFi 



-H-r +1 'iFr 
-H-r+2 iFR-i 



iF 3 
iF 2 
iFi 
H R - X 

H 2 
Hi 



1 



cos k — h 



(cos k — h) 2 + 7 2 sin 2 k 



It is useful to note that Fi-j 
tions, the Toeplitz matrix 



-Fj-i and H t 



Mi = 



-iF H 
H T +iF 



l-j = Hj-i. According to usual conven- 
= Mi[<j>] , M 2 = Mi +M 



in equation (f34t is generated by matrix-valued symbol (analytically continued to the 
unit circle) 

HZ) = ( h{z { -] if(z% 
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while M 2 = Mi - M with 
/ 



M 



-Hi 



V H- R + iF R ■■■ H-i + LFi ff«_i + 
Now, since the trace norm of Mg is vanishing, when R 








. . H-! + LFi 
oo we get dot Mi 



- LFi 

-Fi - iFi 




dct M 2 so that so that C X {R) ~ — 2\/detMi. In order to compute the determinant, 
we need to check the possible zeroes of 

i j., s 1 I i / \ z 2 + l-2/iz 

det c^z) = - 1 + sign(z)^= =_ . , 

1 \ ^(z 2 + 1 - 2hz) 2 - 7 2 (z 2 - l) 2 / 

for instance, when z = ±1, det</>(±l) = | (1 + sign(±l)sign (1 =p /i)) so that for 
/i < 1 the symbol is singular at z = —1, while for h > 1 it is singular at z = +1. 
Unfortunately, as discussed also in ref. l27l , known results for matrix- valued symbol 
do not cover the case of singular symbols with vanishing determinant. Hence, the 
strategy is to factorize the determinant of Mi as a product of determinants of matrices 
generated by scalar-valued symbols. Fortunately, in this case this task is accomplished 
by transforming Mi through the matrix 



U 



iF^H 
1 



so that 



U T MiU 



-iF 

iF 



iH T F !H 



(35) 



Now, we first use a theorem by Widom and Silbermann (see, for instance, 1371 [38)) 
according to which F _1 is a Toeplitz matrix generated by J -1 (in the present case this 
result holds for even i?) and then express the product H T F _1 H as another Toeplitz 
matrix generated by the symbol ih(z~ 1 )f~ 1 (z)h(z). The last identification can be 
done by using repeatedly a theorem by Brown and Halmos 



T((p)T(tp) is a Toeplitz operator iff either tp* (z) or ip(z) are analytic 
functions; if the latter condition is satisfied then T(ip)T(ip) = T(ip?p) 

where T(ip) denotes the Toeplitz matrix generated by the function ip. 
Let us start by computing det(— iF): 

• Haldane phase h < 1 (A 2 > 1). From the analytic continuation to the unit circle 

7 . , , (l-z)(l + z) 



/(*) 



2(1-7) 



sign(z) 



y/(z-X 1 )(z-X 2 )(l-zX 1 )(l-z\ 2 ) 



(36) 
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we see that / vanishes at z = ±1 and is singular at z = A^J,- This case is covered 
by the Fisher-Hartwig conjecture (see refs. l37l [38l l23l E7lD and the asymptotic 
behaviour turns out to be det(— iF) ~ E^R -1 exp(— /3nR) with 

l3 H = 7 -H\X 1 \X 2 ]-ln^— y (37) 

In general the exponent of the power- law prefactor is given by a r ~ Pr where 
the index r runs over all zero and singular points z r of f(z) while the numbers 
a r and j3 r are defined through the factorization of the function in the following 
form 

-vw-rwn(i-i) (i-r 

the residual function t(z) satisfying the conditions of Szego's theorem. In this 
specific case, the 1/R prefactor is due to the combination of the exponents re- 
ported in table [5] Finally, the constant prefactors E can also be written down 
explicitly in the framework of the Fisher-Hartwig conjecture, although the result 
will not be given here for the sake of brevity and because they will be treated as 
free fitting parameters. 

• Critical line h = 1 (A2 = 1)- There is only one zero at z = — 1 and one singu- 
larity at z = 1/Ai. The exponents associated with these two points are the same 
as in the Haldane phase; in this case the power of R receives contributions only 
from these two points and becomes [(1/2) 2 - (3/4) 2 + (-1/4) 2 - (-1/2) 2 ] = 
— 1/2, instead of — 1 . However, the characteristic inverse scale in the exponential 
is nonvanishing even at the critical point 

^ c = ^ln|A 1 |-ln^ Ty . (38) 

• Neel phase h > 1 (A2 < 1). The zeroes remain at z = ±1 while the singularities 
now are at z = A^ 1 and z = A2. Therefore, we proceed along the same line 
followed for the Haldane phase, just by replacing A2 *-* I/A2. In particular, we 
find the same numbers a r and (3 r as for the case h < 1 and thus the asymptotic 
behaviour remains of the form det(— iF) ~ E^R^ 1 cxp(— (3^R) with 

Let us now pass to dct G, with G = iF + iH T F _1 H, generated by the symbol 

1 1 



g(z) = ih(z- 1 )f- 1 (z)h(z)+if(z) 



7 z 2 — 1 

(analytically continued to the unit circle). 



z 2 



2hz + 1 + z^Jiz + z- 1 -2hf 
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Table 5: Values of z r , a r , (3 r for the function f(z) (eq. d36l )) in the Haldane and Neel 
phases. 







0r 


-1 


1/2 


3/4 


+1 


1/2 


3/4 


Ar 1 


-1/4 


-1/2 




-1/4 


-1/2 



• Haldane phase h < 1 (A2 > 1). The Fisher-Hartwig conjecture now can be 
applied, thanks to the following factorization 



g(z) = r{z) (1 - z) 
with 2 and ft^ as in tableland with 



+ ft (1 - (1 + z )™+& (1 + z -i\^-^ 



(40) 



t(z) 



1 



7(i+«r 



- 2/iz + 1 + zJ (z + z- 1 - 2/if - 7 2 (z - z- 1 )^ 



satisfying Szego's theorem. Consequently, the asymptotic behaviour is purely 
exponential: det G ~ E^ exp(— (3#R) where 



1 
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dfcln 



n 



cosfc — h + \l (cosfc — h) 2 + (7sinfc) 2 



+ lnl (41) 



• Critical line h = 1 (A2 = 1). There are no singularities and a simple zero 
at z = —1, with exponents a and j3 as in the first row of table [6] Therefore, 
the net power of R in the algebraic prefactor vanishes and the decay is purely 
exponential with 



1 

2W 



dfc In 



; k — 1 + \/ (cosfc — l) 2 + (7sinfc)' 



(42) 



As a function of 7, /3£ is decreasing for 7 > 1 but does not vanish. 

• Neel phase h > 1 (A2 < 1). With respect to the Haldane phase, the function 
t(z) changes to 



r(z) 



1 



1 



7(l + z) 2 (l-z) 2 



z A - 2hz + 1 + zJ (z + z- 1 - 2/i) - 7 2 (z - z" 1 )' 



while the exponents 011,2 and /3i 2 are reported in the fourth and fifth column of 
table|6] Again, there is no algebraic prefactor and the constant of the exponential 
decay, det G ~ E^ exp(—(3^R), reads 



A - -5F 



2ji 



dk In 



cos k — h + \l (cos k — h) 2 + (7 sin fc) 2 



+ln^=^(7,/i). 
(43) 
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Table 6: Values of z r , a r , f3 r for the function g{z) in the Haldane and Neel phases. 



z r a r (H) /3 r (H) 


a r (N) 


/?r(N) 


-1 1/2 1/2 


1/2 


1/2 


+ 1 -1/2 -1/2 


1/2 


1/2 
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